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An effective field tlieory calculation of the np d^y cross section accurate 
to 1% for center of mass energies 1 MeV is presented. At these energies, 
CN ' which are relevant for big-bang nucleosynthesis, isovector magnetic transitions 

Mly and isovector electric transitions Ely give the dominant contributions. 
The Mly amplitude is calculated up to next-to-next-to-leading order (N^LO), 
and the contribution from the associated four-nucleon-one-photon operator 
is determined from the cold neutron capture rate. The Ely amplitude is 
0^ ■ calculated up to N'^LO. The four-nucleon-one-photon operator contribution 

to Ely is determined from the related deuteron photodisintegration reaction 
'yd np. 
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I. INTRODUCTION 



Big-bang nucleosynthesis (BBN) is a cornerstone of big-bang cosmology [|^. Primordial 
deuterium production is very sensitive to the baryon density of the universe and thus the 
BBN prediction of deuterium abundance can be used to infer this baryon density. These 
deuterium abundance calculations use the cross section for np d'-y as one of the inputs. 
Thus, an accurate estimation of the np — > d'y cross section is essential to the BBN predic- 
tion of deuterium and other light element abundances. At the energies relevant for BBN, 
0.02< E< 0.2 MeV, this reaction is not well-measured experimentally and there are signif- 
icant theoretical uncertainties For example, in the model calculation of Smith, Kawano 
and Malaney (SKM) an error of 5% was assigned to the cross section for np d'-f 0. The 
SKM result agrees with a slightly earlier evaluation from ENDF/B-VI 0. The errors in 
this calculation are not well documented and they could be as large as 10-15% 0. There 
is a much earlier calculation (1967) by Fowler, Caughlan and Zimmerman |^ which agrees 
with the SKM and ENDF/B-VI result for energies 0.1 MeV but disagrees significantly 
at higher energies [0. A recent model independent calculation by Chen and Savage, using 
a low energy effective field theory (EFT) predicted a theoretical uncertainty of 4% . This 
work is a one higher order calculation in the perturbative expansion of the EFT used in 
Ref. 0. The theoretical uncertainty is estimated to be < 1% for center of mass energies 
E< 1 MeV. 

For thermal neutrons {E ~ 10~^ MeV), the cross section for np d'~f is dominated by 
the isovector magnetic transition Mly from the ^5*0 isovector channel to the ^5*1 isoscalar 
channel. At higher energies, E 1 MeV, the cross section for np — > dj is dominated by 
isovector electric transitions Ely from the isovector P-wave to ^5*1 channel. At the energies 
relevant for BBN, the Mly and Ely transitions give comparable contributions. 

Effective field theory is a useful tool in the study of physical processes with a clear 
separation of scales. This is the case for the reaction np — > d'y. At energies relevant for 
BBN, the nucleon center of mass momentum p ~ 30 MeV, the deuteron binding momentum 
7t ~ 45 MeV and the inverse of the singlet channel ^5*0 neutron-proton scattering length 
1/a^^^ ~ 8 MeV are all much smaller than the mass of the lightest meson, the pion with 
mass ~ 140 MeV. Thus a low energy EFT can be constructed for momenta much below 
the pion mass mT^, by integrating out the pions and other heavier degrees of freedom from 
the theory. The strong interaction of the nucleons is then described by four-nucleon local 
operators . The effects of the particles that were integrated out of the theory are encoded 
in a perturbative expansion of local operators, where the expansion parameter is expected 
to be Q/m^r with p, jt, ~ Q. The perturbative description of the low energy physics 

then allows a systematic estimation of errors at any order in the perturbation. 

Recently there has been much discussion about the role EFT can play in im- 

proving the 'traditional' results obtained using effective-range theory (ERT) for low energy 
observables in the two nucleon system. It was shown in Ref. [1^ that for nucleon-nucleon 
scattering, ERT and low energy EFT are equivalent. In some processes, however, involving 
external currents, e.g. electron-deuteron scattering, the ERT amplitude differs from the EFT 
result due to the absence of two-body current operators. The most general set of allowed 
multi-nucleon-external-field operators contains operators that need not be related to the 
nucleon-nucleon scattering operators, included in ERT, by gauging the derivatives through 
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minimal photon coupling. Including multi-nucleon-external-field operators, e.g. two-body 
currents, is straightforward in EFT. The deviation of the ERT result from EFT due to the 
absence of certain two-body current operators will be greater if these operators appear at 
lower order in the perturbative EFT expansion. For example, in the EFT without dynamical 
pions a four-nucleon-one-photon operator contributes to deuteron quadrupole moment iiq 
at next-to-leading order (NLO) [§. The absence of such an operator in the ERT could be 
responsible for potential models' under prediction of /ig by about 5% [§|,|TT| in the impulse 
approximation (lA) ||12|,|13[. In np —>■ dj, the Mly transition amplitude also involves a 
four-nucleon-one-photon operator at NLO in the EFT that is not included in the ERT. This 
operator contributes about 5% to the cold neutron capture cross section a'^^^ = 334.2 mb 
and about 2% at energy E ~ 0.5 MeV [|. On the other hand, the Ely amplitude can be 
written entirely in terms of nucleon-nucleon scattering operators up to N^LO and reduces 
to the ERT result at this order. 

The Mly transition amplitude has been calculated up to NLO PJT^,ITB[. The unknown 



coupling Lf^^, associated with a four-nucleon-one-photon operator, appearing at this order 
can be determined from the cold neutron capture rate of a^^'^ = 334.2 ±0.5 mb [|I^ at 
incident neutron speed v = 2200 m/s. The Ely transition amplitude has been calculated 
up to N^LO 10. In this work, we calculate the Mly transition amplitude up to N^LO and 
the Ely transition amplitude up to N^LO. For the Mly amplitude, there is a new unknown 
coupling L^^^^ at N^LO. Only a linear combination of Lf^^^ and L^^^^ contributes at very 
low momentum. We derive perturbative constrains |ll6| , p!8[| on these couplings to reproduce 
the low energy cross section cr^^^. For the Ely transition, there is a coupling Lf at N^LO 
that is not determined from nucleon-nucleon scattering data. We determine Lf^ from a 
fit to data |jl9| for the related process of deuteron photodisintegration ■yd np. 

The organization of the paper is as follows: We first describe the relevant Lagrangian 
and the power counting rules in Section |I[ This section is rather technical and primarily 
used to define various parameters that enter the expression for the cross section. The 
calculation of the total cross section is presented in Subsections [III A| , [III B| and ^11 C[ The 
renormalization group (RG) flow of the couplings is discussed in some detail, and from 
the analysis, constraints on some of the couplings are derived. In Subsection [III D[ , all the 
remaining couplings are determined. In Section we tabulate the calculated cross section 
for some energies relevant for BBN, discuss the theoretical errors, and compare our results 
with the corresponding values from the on-line ENDF/B-VI database [Q. Summary and 
conclusions follow in Section |V[ 



^We would like to mention that for thermal neutrons, meson exchange currents explain the 10% 
discrepancy in the cross section for np d'j between potential models in lA and a^^^ |14|. How- 
ever, modern potential models, e.g. Argonne vis, with meson exchange currents and relativistic 
corrections underpredict by about 4% |13]. 
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II. THE LAGRANGIAN AND POWER COUNTING 



A perturbative EFT calculation requires a Lagrangian, and a set of power counting rules 
that determine the relative sizes of diagrams contributing to a physical process. The nucleon 
kinetic energy term is 



iDo + 



^0 



2M. 



N 



2M. 



N 



where the covariant derivative is 



1 + 



N 



(2.1) 



(2.2) 



N is an isodoublet field representing the nucleons, the matrix acts on the nucleon isospin 
space and = 938.92 MeV is the isospin averaged value of the nucleon mass. Note 
that the Dq term in Eq. ( |2.1| ) includes all the relativistic corrections to the nucleon kinetic 
energy [g. 

The two-body Lagrangian relevant for radiative capture process np — > d'j and the photo- 
disintegration of the deuteron jd —>■ np can be divided into (a) a Lagrangian £f contributing 
to nucleon- nucleon scattering in the ^Si and ^5*0 channels, (b) a Lagrangian £f contributing 
to nucleon-nucleon scattering in relative angular momentum states ^Pq, ^Pi, ^P2 and (c) a 
Lagrangian /^f^^ describing two-body currents that are not contained in the previous two 
Lagrangians. 



A. Lagrangian for 5'-State Nucleon-Nucleon Interaction: /Zf 



We start by describing the Lagrangian contributing to the final or initial state nucleon- 
nucleon interaction responsible for binding the deuteron. In the '^Si channel, up to N^'LO: 
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-Cq{N^ PiN)\N^ PiN) + y [{N^Of^N)\N^PiN) + h.c. 
-C^{N^OfN)\N^OfN) - ^ [{N^OfN)\N^P,N) + h.c. 
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{N^Ol'^NYiN^Ol'^N) + h.c] + ^ [(iV^Of iV)t(iV^P,iV) + h.c] (2.3) 



-C8°(iV^C»f iV)t(iV^Of)iV) 



c 



{N"'0fN)\N"'0f^N) + h.c 



where summation over the repeated indices is implied. The Pi matrices are used to project 
onto the ^5*1 state. 
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0-2 CTi (&> T2 



Tr[P,tp. 
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(2.4) 



where the a matrices act on the nucleon spin space and the r matrices act on the nucleon 
isospin space. The Galilean invariant covariant derivative operators are defined as: 
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+15 D PiD -6DPiD +PiD 



2 ^ — » , 2' 

D Pi - 2 D Pi D +Pi D 



D Vi - 4 Pi D +6 Pi -4 D Pi +Pi 

■^6 ^5 ^ ^4 ^2 ^3 ^3 

D P^-Q-D PiD+15D P^D -20 D P^T> 



(2.5) 



D D P - D P D - D P D +P D D 



n — 1 



.^ij^xy 



where n is the number of space-time dimensions. For the Mly transition there is also initial 
or final state nucleon-nucleon interaction in the ^5*0 channel and the Lagrangian has the 
same form as the one described in Eq. (|2.3|) , with the corresponding projections onto the 
^5*0 channel. Note that there is no corresponding S-D mixing term, C2'^\ in the ^5*0 channel. 

In 1997, Kaplan, Savage and Wise (KSW) formulated a systematic power counting for 
calculating the cross sections of two-nucleon processes in the context of EFT . One key 
ingredient in the KSW scheme is the non-perturbative treatment of the large ^-channel 
nucleon-nucleon scattering length. KSW power counting can be implemented using power 
divergence subtraction (PDS) where one subtracts the loop-divergences in both n = 3 and 
n = 4 space-time dimensions. An alternative to PDS is off-shell subtraction, where the 



divergences are regulated by subtracting the amplitude at an off-shell momentum point pO 
In this paper, we will use PDS scheme. 

The KSW power counting is as follows: The expansion parameter is Q/A. The nucleon 
center of mass momentum p, the deuteron binding momentum jt, the inverse of the ^5*0 
channel scattering length 1/a^^'^ and the renormalization scale fi are formally considered 
0{Q) and A ~ m^/2 for this low energy EFT (the factor of half comes from the analytic 
structure of the one-pion-exchange contributions). The couplings C2n scale as l/Q"^^, C2n 
scale as and 02^^^ scales as 1/Q. For the low energy theory, we formally take rriT^/Mj^ ~ 
Q/A. Thus, relativistic corrections which come in as p^/M^, ^IjM'^ ~ Q"^ /Mfj = Q"^ /m\ x 
ml/M% contribute at N^LO. 

It is a feature of the KSW power counting that the EFT couplings associated with S- 
state interactions scale with some powers of Q. This is because the couplings are fine-tuned 
to reproduce the large scattering lengths that one sees in the ^5*0 and ^5*1 channels. In order 
to reproduce the exact deuteron pole, the coefficients C2„ are expanded in powers of Q: 



Co 

C2 



Co,-i + C'o.o + ■ ■ 
C2-2 + C2-1 + ■ 



(2.6) 



where the second subscript denotes the Q scaling. The EFT couplings C2n are determined 
by matching the nucleon-nucleon scattering amplitude calculated in EFT 0,^,|1^] and that 
obtained from the ERE : 
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pcotS = -7t + {f + -it) + ]^W2 (f + I'ff + ■ 
To the order we are working, the deuteron binding momentum 7^ 



7 = y/M]\jB, and B = 2.2246 MeV [E^ is the deuteron binding energy. The couphng C2 
3r]d/ {V^l"^) is determined from the S-D mixing parameter ei p|, p!6| , p3| , where rjd 



(2.7) 

7 - 7V (8M^), where 

(sd) _ 

2 — 
0.02534 



is the D-wave to 5*- wave ratio at the deuteron pole p2|]. Notice that €2^^^ is smaller than 



naive power counting estimate because of the smallness of the parameter r]d- 

For our calculation, we need the ^5*1 channel initial or final state interaction up to N^LO 
and the ^Sq interaction up to N^LO. In the ^Si channel, at N^LO only four experimental 
inputs (7, Pd, W2 and rjd) enter the scattering amplitude. The EFT couplings that depend on 
the high energy scale can all be expressed in terms of these four parameters. However, only 
the combinations C4 + C4 and Cg + Cg = Cg contribute to nucleon-nucleon scattering and 
so nucleon-nucleon scattering data cannot be used to separate C4 from C4, and Cg from Cg . 
We use the effective range pd = 1.764 fm and shape parameter W2 = 0.778 fm^ |^^. In the 
^5*0 channel, at N^LO the nucleon-nucleon scattering amplitude is completely determined by 
the scattering length a^^^ and the effective range tq. The EFT couplings are determined 
in terms of these two parameters. Note that the experimentally measured scattering length 
Q^exp^ for neutron-proton scattering in the singlet channel ^5*0, gets about a 3% contribution 

0. 



from magnetic moment interaction 



Isospin-breaking single magnetic photon exchange 



gives the dominant correction and it is described by the Lagrangian 



2M" 



N 



(2.8) 



where B is the magnetic field and ki = 2.353 is the isovector nucleon magnetic moment 
in nuclear magnetons. These N^LO interactions must be included in the Mly amplitude. 
We implicitly include the magnetic moment interactions by using the experimental value 



exp 



—23.75 fm For the effective range we take tq = 2.73 fm. 



B. Lagrangian for P-State Nucleon-Nucleon Interaction: £2 

In this subsection, the second Lagrangian £^ is described. For the Ely transition 
amplitude, there are also initial or final state nucleon-nucleon interactions in relative angular 
momentum states ^Pq, ^Pi,^P2 described by the Lagrangian [H: 



'~'2 — I ^2 "xy"wz \ ^2 \Pxw(jyz "xz'Jyw] T ^2 



4 

'^^xw^yz ~l~ '^^xz^yw ~^xy^wz 



xl(iV^O(f )iV)t(iV^Oa/)iV) , (2.9) 
where the P-wave operator is 

dif' = h.lf'-lf'^., (2.10) 
and Pj are the spin-isospin projectors for the isotriplet, spintriplet channel 
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Pi = —r^cr2Cri T2T3 



Tr 



p(p)pipy 



—5j^ ■ 
2 ^ 



(2.11) 



The power counting of the P-state couphngs is straightforward. Since there are no fine 
tuned high energy scales in the P-state scattering, the couphngs that are dependent only 
on the high energy physics are 0{1) i.e. they do not scale with Q. As a result, the next set 
of P-state operators, suppressed by two extra powers of momentum p, enters only at N^LO. 
Only the linear combination 
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(2.12) 



enters our calculation and Nijmegen phase shift analysis |T2| , p5| fixes Cp = —1.49 fm*^ 



C. Lagrangian for Two-Body Currents: £2^'^ 



Finally, we discuss the two-body currents contributing to np — » d'y or 7c/ np, which are 
not included in the previous two Lagrangians £f and These operators are not related 
by gauge transformation to the nucleon-nucleon scattering operators in £f and /^f. 

-eZf ^^(iV^P,iV)t(iV^0f iV)B, + ^eLf ^(iV^o£'^)iV)t(iV^P,iV)E, 

-leLf -(iV^0£'^)iV)t(iV^0f iV)E, - L^'^ {N^ Ol'f^ N)^ {N^ PaN)E, (2.13) 

+ ^eLf ^(iV^0£'^)iV)t(iV^0(4)iV)Ei + ^eZf ^(iV^Og'^)iV)t(iV^Op)iV)E, 
-eL^'Hei,kiN''PiNy{N^PjN)Bk + h.c. , 

where E and B are the electric and magnetic fields respectively, and 

1 



Pi = —;=cr2 ® T2Ti 



(3,P) 



■^2 

D P. 



2 D Pi D +Pi D 



(2.14) 



DiD Pf ^ - 2 D.Dfc Pf ^ Dfc 



(P) 



D. PP D 

D PP D. +2 hk PP DfcD. -PP d'd. 



(P) 



The superscript on the two-body current couplings L's denote the transitions that the par- 
ticular operators contribute to. For the Mly transition only a specific p independent com- 
bination of Li^^^ and L^^^^ contributes and we fix it from the cold neutron capture rate 
np d'y. Renormalization group (RG) flow analysis determines Z^^^^ in terms Lf^^^, C2 "^"^ 
and C2 ■ For the Ely transition, only a combination, L^i, of Lf^^, Lf^^ and L^^^ enters 
the calculation at N^LO and it is fixed from a fit to the related deuteron breakup process 
'jd —>■ np The other operators contribute at orders higher than considered here. The 
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RG analysis of these operators is discussed in more detail below, in Subsections |III A|, \Lll B 



Finally, the observed value of the deuteron magnetic moment hm fixes the isoscalar magnetic 



'(Mis) 



-0.149 fm^ g,^ 



moment coupling L2 

In the last three subsections, we described the Lagrangian and the power counting rules 
relevant to our calculation. Now, the calculation for the total cross section is presented, along 
with the estimated theoretical uncertainty. This is followed by the discussion of a matching 



procedure for determining the unknown couplings 
fixed from the jd —>■ np data. 



MU 



and L 



Mil 



The parameter Lei is 



III. CALCULATION AND ANALYSIS 



The amplitude T for low energy np — > ^7 reaction is [|1^,|TB]: 



+e Xmis-^U T2 0-2 

1 J. 
+e Xe2s-^U T2 0-2 



u 



(3.1) 



U , 



where only the lowest partial waves are shown: isovector electric dipole capture of nucleons 
in a P-wave with amplitude XEiy, isovector magnetic dipole capture of nucleons in ^Sq 
state with amplitude Xmiv, isoscalar magnetic dipole capture of nucleons in ^Si state with 
amplitude Xmis and isoscalar electric quadrupole capture of nucleons with amplitude Xe2s- 
We use dimensional regularization to regulate divergences and n represents the number of 
space-time dimensions. The t/'s are two component nucleon spinor wave functions. \p\ =p\s 
the nucleon center of mass momentum, k is the outgoing photon momentum, and e(^) and e(d) 
are the photon and deuteron polarization vectors respectively. The following dimensionless 
amplitudes, X, are defined 



\p\M, 



N 



7^ + p^ 



X 



EU 



Xmic 



2 



Xeu 



—Xmi 



Xmu 



2 



Xmu 



and the total cross section for np — * d'j is then written as: 
^ 2/ + ^ 3^4 



—Xe2s ) 



(3.2) 



4M2,(p2 + ^2) J ^^M%p 2 



dcosO 



3, ~ 
2 



XeivI 



sin^ 9 



\^MU 



+ \Xm1s\'^ + I^S2sl 



(3.3) 



2p^ + 4pV + 37^A (7^+p' 



2^3 



{XeIv^ + i^Mlv)'^ + (-^Mls)^ + {^E2s)^ 



\ 4M2,(p2+^2) J ^^M%p 

where the deuteron mass is ^ 2Mjv — jM^q The relativistic corrections enter 
at N^LO and so only the E\y cross section which is calculated to this order receives such 
contributions. 
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A. Isovector Electric Transition: Ely 



The Ely amplitude up to N^LO is 



MnPI^ ^ Pd7 ^ 3 2 2 5 3 3 35 4 4 5 2 , (p' + 7')% 



[pi + J 



2\2 



+—75- cos y — 



16M2,(p2 + ^2 

where the relativistic corrections contribute at N^LO as they are suppressed by additional 
powers of (m.^/Miy)'^, and 



'0,-1 



Lei = -(35^) 13. 5j 



■^'-^4,-1 = '-^4-1 "T 



(/i - 7)3 



The TpA,-! coupling renormalizes contributions proportional to 77^ in Eq. ( p.4|) . The renor- 
malization scale /i independent parameter L^i ~ can be simplified further. We assume 
that all the couplings have a sensible large scattering length (7 0) limit. Then it follows 
that at the high scale /z = A, only the combination (LfLi — MN-rpi,-i)/C^ ^{^ contributes 

a 0(1) term. The other combinations, e.g. 'j'^Lfl^^/CQ contribute terms of 0{Q^) and 
higher. Thus RG analysis gives 



-^1,-1 ~ M^rpi-l 
'0,-1 



Lei - ' ^(3,;;" "'^ ( 1 + 0{Q'^) ) , (3.6) 



and we keep only the leading contribution to Lei for our calculation. It is possible to make 
an order of magnitude estimate for Lei- The contribution from S-D mixing, proportional 
to ?7^, is numerically negligible at N^LO even though it is formally a N^LO term. Thus, 
ignoring contribution from ^C4^_i, we estimate at the renormalization scale = A ~ 

47r 1 . 1 

l*-^o-i I ~ H~T' 1-^1,-1 



MnA MnA^ 
=^ |Lsi|~2fm^ (3.7) 

where factors of 1/Mn from the nucleon loops were included. 

In Eq. ( p.4| ), we have only kept the terms that contribute to the total cross section after 
summing over deuteron, photon and nucleon polarizations. We have included the formally 
N^LO contribution proportional to cos^. This term, odd in cos^, contributes to the cross 
section only at N^LO after the integration over the angle 6, Eq. ( p.3|) . In obtaining Eq. ( p.4|) , 
we have also used the fact that only the fj, independent combination 
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enters our result. From RG analysis, by flowing the couplings to the high scale /x = A, 
we see that the numerator is 0{1/Q) instead of the naive counting of Thus this 

particular combination of operators contributes starting at N^LO and we ignore it in this 
N^LO calculation. 

The dominant contribution beyond leading order, Eq. (|3.4|), is a simple expansion of the 



(3.9) 



factor A/Zrf = 1/ \/l — pdl from the deuteron wave function renormalization 

8717 



Z 



1 , ,22,33,44, 2 ^7 - 5/i 7^ 7/i - 57 
I + Pdl + Pdl +Pdl +Pdl +Vd— — + ^7J2— — + 



p-'y SMf, /i - 7 



In Eq. (|3.9|) , the pdl terms inside the square brackets are from the expansion of Zd-, and 
the sixth and seventh terms arise from the S-D mixing and relativistic corrections to the 
deuteron two-point function respectively. 

Amplitudes for inelastic processes, e.g. Ely and Mly transitions, with a deuteron in 
the final state involve a factor of a/Z^. On the other hand amplitudes for elastic processes 
involving a deuteron, e.g. the deuteron quadrupole form factor Fq, include a factor of Zd- 
In Ref . , an expansion in Zd — 1 was proposed which resums the expansion in pdi seen in 
Eq. ( |3.9| ) and thereby reproduces the exact factor of Zd at NLO for the amplitude for elastic 
processes. On the other hand, amplitudes involving inelastic processes on the deuteron 
have this sum incomplete. The factor of y/Zd is reproduced only in perturbation. This 
incomplete sum can be avoided if we do an expansion in a/Z^ — 1 to reproduce the exact 
factor of y/Zd for inelastic processes. The overall factor of Zd in elastic processes is then 
reproduced at N^LO instead of NLO. Thus, one might think that the a/Z^— 1 expansion is 
more appropriate to calculate the cross section for the inelastic process np — >■ di. However, 
for the total cross section, which is summed over the nucleon spins and the deuteron and 
photon polarizations, it is the square of the amplitudes that enter the expression. In such 
cases it is more convenient to do the expansion in — 1 and reproduce the exact factor of 
Zd at NLO. The three different expansions pdi, Zd — 1, v^^— 1 are formally equivalent but 
correspond to different ways of relating the C2n couplings to the ERE. 

In the Zd — 1 expansion, the position of the deuteron pole at momentum p = ijt ^ il 
and the residue Zd (without S-D mixing) at the pole are used to fix the C2„ couplings. One 
can think of the Zd — 1 expansion as reproducing the asymptotic tail of the deuteron wave 
function exactly. This allows a better description of data as low energy experiments probe 
only the large distance behavior of the deuteron wave function. For the rest of the paper we 
will use the this expansion. The results for X can then be written as: 



xei.= - /r:"'o.. 1 + ^{Zd - 1) - ^{Zd - ly + i^iZd - If - —{Zd - ir 
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ip' + ir 



127r \ 3 / V 2' V 2 2M, 



-L 



El 



N 



8^4 + 17^2 2 + 54 1^^ 1 ,.2\^P' 2n 

-iZd- I) ]+ ^TTT COS 9 
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{XeivY 




(3.10) 



where the last part in {XeivY is from the cos^ 9 terms. 

The expression for (Xeiy)"^ in Eq. (|3.10| ) is a perturbative resuh. There are corrections 
to this expression from terms that are higher order than the ones considered here in the 
perturbative expansion. The higher order corrections can be separated into: (a) factors of 

from the wavefunction renormahzation of the deuteron, Eq. and (b) contributions 
from higher order nucleon-nucleon scattering operators and two-body currents. At ~ 1 
MeV, we expect higher order corrections to contribute {Z^ — l)(p^ + 'y'^)/Mjf ~ 0.004, 
5{Zd — 1)?7^ ~ 0.002 due to wave function renormahzation. The Z^ — 1 corrections to 
the term involving Lei are renormalized away when we fit Lei- The higher order P-wave 
operators should contribute 7p^/(67rA^), 'j^p'^ /{6nA^), 7^/(67rA^) ~ 0.01 with A ~ m,T^/2. 
Thus we estimate the theoretical error to be ~ 1% at energies oi E 1 MeV. The error is 
smaller at lower energies. 

B. Isovector Magnetic Transition: Mly 

A straight forward calculation of Mly amplitude gives: 



X 



(0) 
MU 



X 



(1) 

MU 



+ 



■'np 



X 



(2) 
MU 



^exp 



IP 



-i(Z,-l)2xlt + ^3-^I 
where the /x independent parameters are 



np } 



(3.11) 



Lnp — 


{fi 


-7)(/i- 




tMIv 
-^1,-2 


Lnp — 


(/^ 


-7)(/i- 


l/a"^P) 


tMIv 
-^1,-1 




(/i 


-7)(/i- 


l/a"^f) 


tMIv 
^1-1 



+ 



M^^\{^i-^Y (/i-7)^ 

2tMIv , '^l^ i^d - 1)^ 
7 ^3,-3 + -T7 7 

MnI (/i - 7) . 

KlTT (Zd - 1 



(l + 0{Q' 



(3.12) 



An analysis similar to the one used for Lei-, allows us to simplify the parameters LnpS above 
in Eq. ( |3.12| ). A naive estimate gives 



Lnp ~ 



-9 fm^ 



Lnp ~ +3 fm^. 



(3.13) 
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We use the experimentally measured scattering length a^^^ which incudes the effect of single 
magnetic photon exchange in neutron-proton scattering. In Eq. (|3.11| ), from RG analysis, 
we used: 



'-Mlv 
^3,-3 



1 



^1.-2 + 



1 



(3.14) 



The first part 0{1/Q'^) gives higher order corrections to 1^^^^ ~ So we can set ^ = 



at N^LO, and have 1^^^^ completely determined. Up to N^LO 



I Mlv 



(3.15) 



The singlet channel scattering length a^^'^ = —23.75 fm is unnaturally large and it is easier 
to do the error analysis for the cross section in the limit a^^^ oo, with finite corrections 
from l/(a*^^^7) and l/{a'^'^^p) terms. There are corrections from initial state interactions in 
the ^So channel that come in as TqP^/'j^ ~ 0.03 at center of mass energy E = 1 MeV. There 
is also a contribution from ^Sq operators that come in as ~ —0.002, where the shape 

parameter ri ~ — 1 fm^ is the ^5*0 equivalent of W2 in Eq. ( |2.7| ). The momentum independent 
corrections from factors of Z^ — 1 from deuteron wavefunction renormalization, ro7 from ^Sq 
channel initial state interaction etc. get renormalized away when we fit the parameters L„p's 
at very low momentum. We estimate the errors in Mly to be ~ 3% aX E = 1 MeV. 



C. Isoscalar Magnetic and Electric Transitions: Mis, E2<. 



The Mis and E2s amplitude have been calculated for cold neutron capture |T5|,p!6 
leading contribution from these transitions at non-zero momentum transfer is: 



The 



\X 



Mlc 



Mlc 



\Xe2c.? 



27r2(p2 + ^2 

iWd {p^ + 7^272 + ^4) 



) 



100 



(p2 _j_ ^ 



2^4 



(3.16) 



The coupling L2 



can be fitted to the deuteron magnetic moment at the renormalization 

Mis I 



scale /i = and gives L2^'^^{'mn) = —0.149 fm^ HHH- 

In the last three subsections, the np d'y amplitude in terms of the lowest partial waves, 
Mlv, Ely, Mis and E2s has been calculated in terms of three unknown parameters Lnp, 



Lnp and Lei, Eqs. ( |3.10| ), ( |3.11| ) and ( ^.l(j| ). We now move on to describe how the three 



unknown parameters Lnp, Lnp and Lei can be determined from np — > d'y for cold neutrons 
and from photodisintegration of the deuteron ■yd —>■ np. 



D. Determining the Parameters Lnp, Lnp and Lei 

The leading contributions from each partial wave are shown in Fig. |l|. One can see 
that the Mis and £"25 transitions can be ignored for a 1% level calculation in the energy 
range of interest. At center of mass momentum po = 0.003443 MeV the cross section for 
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FIG. 1. np dj cross section a (mb) for center of mass energy E (MeV), on a log-log plot. 
Solid curve: LO Mly transitions, dash: LO Ely transitions, long-dash: LO Mis transitions and 
dot-dash: LO E2s transitions. 



np — » d'y is measured to be a^^^ = 334.2±0.5 mb |Tj] and Mly transitions give the dominant 
contribution. Fig. |T]. Thus at NLO, we fix Lnp from this cold neutron capture rate: 



(LO) _ 47ra(7'+p')^ ^(o) ,2 



a 



(NLO) 



P=PO 
2\3 



4na{-f^+p-r ^„ .y(o) .^(1) 



■'np 



9.039 ± 0.027 fm^ • 



V=PO 



(3.17) 



The N^LO terms also contribute to the cross section at momentum po- Therefore, we impose 
the constraint 



(NNLO) — 



47rQ;(7^ + p 



,2\3 



Ko) 



'(2) 







P=PO 



to reproduce the experimental cross section ct'^^p. This determines the parameter 



-'np 



4.957 ±0.011 fm' 



(3.18) 



(3.19) 



The parameter Lei, which contributes only to Ely transitions, cannot be determined re- 
liably from this single data point since the Ely cross section is negligible at this low mo- 
mentum pq. The error due to the experimental uncertainty in cr'^^^ is significant for energies 
E^ 0.25 MeV. Above these energies the Ely cross section gives the dominant contribution 
and the error in Mly cross section can be ignored. The contribution of the four-nucleon-one- 
photon operator is found to be significant at low momentum pq = 0.003443 MeV. Neglecting 
this operator by setting L^^^ = at NLO underpredicts the experimental result cx^^'^ by 
5%. This number is different from the one quoted in Ref. [^,15|, where the pdj expansion 
was used. 

The parameter Lei can be determined from photodisintegration of the deuteron. This 
cross section is related to the neutron capture cross section by ||6|JT9[|: 
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2.5 

2 

a(mb) 1-5 
1 

0.5 



7(i — > np 



10 



E^(MeV) 



FIG. 2. Cross section in mb for jd — > np as a function of photon energy in MeV. The data 



is taken from Ref. |19], pages 78 and 79. 



a{'jd np) 



2M, 



N 



3^2 



— B 



6BE^ + AE^ 



4M, 



TV 



a{np d'j) 



(3.20) 



where E^ is the incident photon energy in the deuteron rest frame. Only the Ely transitions, 
on the right hand side of Eq. ( |3.20|) , receive relativistic corrections. In principle it is possible 
to determine all the three parameters Lnp and Lei from the deuteron breakup data. 
However, such a determination leads to significant uncertainty in the Mly cross section. 
This is because there are only two data points in the nucleon energy range sensitive to 
Mlv transitions, 0.3 MeV < E k. E^ ~ B< 0.5 MeV. We find that a fit to data p 



m 



the photon energy range 2.6 MeV < E^ < 7.3 MeV does not give a reliable constraint on 
the parameters Lnp, Lnp and Lei- Constraining Lnp and Lnp from a'^^^ as described above 
determines them more accurately. The other parameter Lei is determined from a fit to 
data in the photon energy range 2.6 MeV < E^ < 7.3 MeV, where experimental errors 
in the total cross section from the Mly transition are negligible. This gives: 



Lei = -(5.3 ±3.6) fm^ 



(3.21) 



The error due to experimental uncertainty is < 1% over the range of E^ 5 MeV. This means 
that for nucleon energies relevant for BBN, E^ 0.2 MeV, the experimental uncertainty is 
significant compared to theoretical error, see Fig. ^ A few more high precision measure- 
ments in the incident photon energy range 2.5 MeV < E^^ 5 MeV would provide important 
constraints on Lei and determine the np d'~f cross section more accurately at energies 
relevant for big-bang nucleosynthesis. 

The contribution from N^LO is found to be ^ 3% for incident photon energies E^'^ 8 
MeV and is a small correction to the N^LO result [0. This is better than the naive theoretical 
estimate for energies E^ > A MeV. This low energy theory, which is formally valid for energies 
E^< 8 MeV, seems to reproduce data well above its range of validity, see Fig. 0. 

Before discussing the results, it is important to check the self-consistency of our pertur- 
bative EFT calculation. The calculations were carried out under certain assumptions about 
the sizes of the couplings. The fact that the fitted couplings, Eqs. ( |3.17|) , ( |3.19|) and ( p. 211) , 
have sizes comparable to the theoretical estimates, Eqs. (|3.13|) and ( p.7|) , suggests that the 
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np d'y 



% Error 
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E (MeV) 

FIG. 3. Estimated theoretical error and experimental uncertainty in percentages versus nu- 
cleon center of mass energy E (MeV) for the np d'y cross section. Solid curve: experimental 
uncertainty, dotted curve: estimated theoretical error. 

perturbative expansion is under control to the order calculated. This allows us to make 
reliable error estimates. 



Table | shows the EFT np — ^7 cross section for various nucleon center of mass energies, 
E. The corresponding values for the cross section from the on-line ENDF/B-VI database p| 
are also shown in the last column. As explained earlier, an error of (rip'^/7 + (1 + 

|l/(a^^^p)| + 1 1/(0*^^^7) I) with respect to LO was assigned to the Mly cross section and 
a conservative error of 7(p^ + p^7^ + 7^)/(67rA^) — h^Zd — l)rjl + {Z^ — + 7^)/M^ 
was assigned to the Ely cross section. The errors were added linearly to the total cross 
section and are found to be ^ 1% for E"^ 1 MeV. The EFT results are presented to only 
four significant digits, unless the theoretical errors enter earlier, in which case we keep up 
to the first digit where the errors contribute. Compared to the results in Ref. 0, these 
higher order EFT results are perturbatively closer to the ENDF/B-VI values. However, at 
some energies, e.g. 1 x 10~^ MeV, where the difference between the N^LO EFT result and 
ENDF value is much larger than the expected perturbative corrections, not surprisingly the 
discrepancy does not disappear by going to N'^LO p. 



In finale, the big-bang nucleosynthesis prediction for primordial light element abundances 
uses the cross section for np — >• d-y as an input. A precise estimate of this cross section is 
thus critical for these predictions. The total cross section for radiative capture of neutrons 
on the proton np ^7 was calculated in EFT for center of mass energies E"^ 1 MeV. At 
the energies relevant for BEN, 0.02 MeV ^ E^ 0.2 MeV, the isovector electric transition 
Ely and the isovector magnetic transition Mly give the dominant contributions. 



IV. RESULTS 



V. SUMMARY AND CONCLUSIONS 
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a{np — > d'y) 



Hj ^^iVlc V j 


iVl 1\/ I lilU I 


J-j lY J 


A/f 1 4- F1 fmVi'l 
IVl 1\/ JZ/ 1^ I lilU I 


rji\Ljr yiiLU J |I0I] 




qq /I 9 * 

0041:.Z 


|^ 1 99 V 1 
O.IZZ X lU 


qq /I 9 * 


ooz.u 


0. X iU 


i.DD < ^Uj 


n nni ni q/'"?'! 
u.uuiUiy(^o ) 


1 f;f;a('n^ 

i.DDo^U j 


1 KRn 

i.DDU 


1. X lU 


1 i7ni'n'\ 

i . i ( Ul^Uj 


n nni /i/ii /'^^^ 




i.iyo 


5. X 10-3 


0.4950(0) 


0.00322(1) 


0.4982(0) 


0.496 


1. X 10-2 


0.3279(0) 


0.00454(2) 


0.3324(0) 


0.324 


5. X 10^2 


0.09810(0) 


0.00997(3) 


0.1081(0) 


0.108 


0.100 


0.04973(0) 


0.01379(5) 


0.06352(5) 


0.0633 


0.500 


0.00787(3) 


0.0263(1) 


0.0341(2) 


0.0345 


1.00 


0.0036(1) 


0.0313(2) 


0.0349(3) 


0.0342 



TABLE I. Cross section for np d'j in mb for different center of mass energy E (MeV). For 
comparison, the values from the ENDF/B-VI on-hne database are also shown. The first entry (*) 
is used for fitting a combination of parameters Lnp's, Eq. (|3.17| ), (|3.19| ). Another parameter Lei 
was fitted to the deuteron photodisintegration "yd — > np cross section, Eq. (|]2l|). 



We calculated the Ely cross section to N^LO and the Mly cross section to N^LO. Up 
to N^LO, the Ely amplitude calculated in EFT is equivalent to the effective range the- 
ory result. Thus in principle, any calculation including potential models that reproduces 
nucleon-nucleon scattering data well will be able to describe the Ely amplitude to the ac- 
curacy of the N^LO calculation. However, the effective range theory result differs from the 
N^LO EFT result due to the absence of a four-nucleon-one-photon operator. This operator 
cannot be obtained from nucleon-nucleon scattering data alone. We determine this oper- 
ator from the deuteron breakup jd —>■ np data with significant experimental uncertainty. 
Similarly, for the Mly transition there is a four-nucleon-one-photon operator at NLO that 
is not related to nucleon-nucleon scattering operators by gauge transformations. We deter- 
mine this unknown operator from the cold neutron capture np — »• rate and this operator 
contributes ~ 2% to the total cross section at E = 0.5 MeV. A few more precise measure- 
ments in the incident photon energy range 2.5 MeV < E^^ 5 MeV would provide important 
constraints on the Mly and Ely transitions in the energies relevant for BBN. The sizes 
of the operators determined from data are consistent with their theoretical estimates. This 
verifies the theoretical assumptions about the perturbative expansion to the order of the 
calculation, which allows us to make reliable error estimates. For energies E< 1 MeV, the 
theoretical uncertainty in the total cross section is estimated to be < 1%. 
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